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Abstract 

We identify a new superintegrable Hamiltonian in 3 degrees of freedom, obtained as a reduction 
of pure Keplerian motion in 6 dimensions. The new Hamiltonian is a generalization of the Keplerian 
one, and has the familiar 1/r potential with three barrier terms preventing the particle crossing 
the principal planes. In 3 degrees of freedom, there are 5 functionally independent integrals of 
motion, and all bound, classical trajectories are closed and strictly periodic. The generalisation 
of the Laplace-Runge-Lenz vector is identified and shown to provide functionally independent 
isolating integrals. They are quartic in the momenta and do not arise from separability of the 
Hamilton-Jacobi equation. A formulation of the system in action-angle variables is presented. 
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I. INTRODUCTION 



The Kepler problem is well known to be superintegrable - that is, it has five functionally 
independent integrals of motion. They are the energy and the components of angular mo- 
mentum and the Laplace- Runge-Lenz vectors, obtainable by separating the Hamilton- Jacobi 
equation in spherical polar and rotational parabolic coordinates [l). For Hamiltonians with 
three degrees of freedom, the existence of five integrals of motion implies that every bound 
trajectory is closed. 

Sommerfeld and Born, in the days of the old quantum theory, appear to have been the 
first to realize that if a potential is separable in more than one coordinate system, it possesses 
additional isolating functionally independent integrals The first systematic inquiry 

into this problem was begun by Winternitz, Smorodinsky and co-workers, who found every 
potential in two degrees of freedom for which the Hamilton- Jacobi equation is separable in 
more than one way [4]. Subsequently, they extended this work to three degrees of freedom 
by finding every potential separable in spherical polars and at least one additional coor- 
dinate system [5|. Evans 6] then completed this work by investigating all the remaining 
possibilities. A useful introduction to the subject of superintegrability, as well as summary 
of recent work, is given in the conference proceedings of Tempesta et al. Q]. 

If the Hamilton- Jacobi equation separates, then the corresponding integral of motion is 
necessarily linear or quadratic in the canonical momenta. Consequently, all the superinte- 
grable systems listed in {4], 0, [f| have integrals that are quadratic in the momenta. As an 
example, let us consider the Keplerian Hamiltonian 

H = -\v\ 2 --, (1) 

where A; is a real positive constants. Letting r denote the position vector, then the integrals 
are the energy E, the components of the angular monetum vector L = r x f and the 
Laplace-Runge-Lenz vector 

k 

A = r x L r, (2) 

r 

all of which are at most quadratic in the velocities. 

Superintegrable systems with higher-order integrals are known jg], although they are 
extremely scarce. Examples include the anisotropic harmonic with rational frequency ratio 
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£ : m : n where £ + m + n > 5, i.e., 

# = |ff + £V + m y + 71V. (3) 

The potential separates in rectangular cartesians and possesses two commuting quadratic 
integrals. There are two additional integrals which may be taken as polynomials of degree 
£ + m — 1 and £ + n — 1 [9] . The Calogero potential in a harmonic well 



1 9 , 9 n 9. A^l /Cl /Cl 

2 [x — yy (y — z) 1 [z — x) 1 



is known to be super-integrable and possesses an integral of the motion that is cubic in the 



velocities 10 



In this paper, we introduce a new superintegrable Hamiltonian, namely 

1i.i9 k ki k 2 k-x . . 

H = | r | 2 -- + 4 + ^ + 4 (5) 

2 r x A y z z A 

This is recognized as a generalization of the familiar Keplerian Hamiltonian. The constants 
hi, k 2 and k 3 are taken as positive, so the Hamiltonian is perfectly physical and motion is 
confined to, say, the octant x > 0, y > and z > 0. The terms involving the ki correspond 
to repulsive barriers preventing the orbit crossing any of the principal planes. In Section [Til 
we present numerical integrations showing that the orbits give closed curves in all cases. 
This motivates a search for the isolating integrals, one of which is found to be quartic in the 
momenta in Section IIHI Finally, an action angle formalism is given in Section IIVI and the 
relation to the Kepler problem discussed. 



II. EVIDENCE OF SUPERINTEGRABILITY 
A. Analytic Proof 

Let us recall that Keplerian motion in iV degrees of freedom always possesses 2N — 1 



functionally independent integrals of motion [12|. Specialising to 6 degrees of freedom, we 
have the Hamiltonian 
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where s has Cartesian coordinates (s±, S2, S3, S4, S5, sq). Now, let us introduce coordinates 
(x, y, z, 8 X , 9 y , 8 Z ) according to 

Si = icos^, s 2 = xsin6' x 

s 3 = y cos 8 y , s 4 = y sin 8 y 

S5 = zcos9 z , sq = z sin 8 Z 

The Hamiltonian becomes 

^ - o tPl + ^2 + P3 + —o + TT + -J2J 



nU'l 1 ^2 1 ^3 1 o „2 1 r 2 ' / 9 1 2 1 2 \ — 

The coordinates (81, 82, #3) are ignorable, so we obtain a new Hamiltonian 



where fci,/c2 and &3 are the constant values of the momenta conjugate to the ignorable 
coordinates. In the original 6 degrees of freedom Hamiltonian (J6j), every bound trajectory 
is closed. Consequently, in the reduced 3 degrees of freedom Hamiltonian (CO), every bound 
trajectory is also closed. Evidently, the proof can be readily generalised to N degrees of 
freedom. 



B. A Sampler of Orbits 

It is interesting to investigate characteristic orbits corresponding to the Hamiltonian (JTj). 
Using a standard Burlisch-Stoer code 13j, the motion of a unit mass particle was followed 
in the octant with x > 0, y > and z > 0. Stepsize and tolerances were set to maintain 
accuracy to a level of around 10~ 12 relative energy change, and integration lengths were 
typically tens of periods. 

Starting with initial conditions that would give a circular orbit in a true Keplerian poten- 
tial, the effect of relatively weak barriers can be seen in Figure HJ The orbit appears similar 
in shape to the Keplerian ellipse but reflects off the three axes planes. In this case, as a 
consequence of the ki being equal and the symmetry in the initial phase space position, the 
orbit is confined to a plane. The effect of larger barriers can be seen through increasing the 
ki by a factor of ten, as shown in Figure [2j The orbit is now further distorted, and takes 
the form of a figure-of-eight. If the ki are not equal, the orbit is still closed, as shown in 
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FIG. 1: A perturbed circular Keplerian orbit in the potential. The three panels show the projec- 
tions in the x — y, x — z and y — z planes respectively from left to right. The initial conditions are 
given (to 2 s.f.) below the plot and the energy in the top corner. Note that although in this case 
the orbit lies in a plane this is not generally true. 




x = (5.8, 5.8, 5.8) p = (0.0, 0.3, 0.0) 



*i = (0.10, 0.10, 0.10) 



k= 1.0 



FIG. 2: As for Figure [T] but now the centrifugal barriers are an order of magnitude larger. 

Figures [3] and 0J The latter of these two cases has a different set of initial conditions, which 
would place it on an initially elliptical orbit in the true Keplerian problem. 

Many more initial conditions and combinations of parameter values were investigated. In 
all bound cases, every orbit is closed and strictly periodic, which is reassuring confirmation 
of the existence of a fifth isolating integral. 
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x = (5.8, 5.8, 5.8) p = (0.0, 0.3, 0.0) k , = (0.55, 0.20, 0.10) k=1.0 



FIG. 3: As for Figure [2] but now the centrifugal barriers are different and the motion is no longer 
confined to a plane. 
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x = (10.0, 10.0, 10.0) p = (0.1, 0.0, -0.2) ^ = (0.10, 0.10, 0.30) k=1.0 



FIG. 4: This orbit in the potential has initial conditions that would place it on an elliptic orbit in 
the original Kepler problem. 

III. THE INTEGRALS OF MOTION 

Makarov et al. jlj] and Evans [(J already showed that all Hamiltonians of the form 

1 . . n fci ko k"i , . 

H = -v 2 + F(r +4 + ^ + 4 8 

2 ar y z z z 

possess four isolating integrals of motion, arising from separability of the Hamilton- Jacobi 
equation in the spherical polar and conical coordinate systems. Here, F(r) is an arbitrary 
function of the spherical polar radius. The four isolating integrals are the energy E and 
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three generalizations of the angular momentum components, namely 

B= -|p|* + F W + -I + -§ + -§ (9) 

/2 = ^ + ^ + ^|! („, 

2 ar 2r 

J 3 = ^3 + ^f + ^ (12) 



h = L lP2 - Pl L 2 -2z (-— + _i + _§) ( 13 ) 



where p and L are the linear momentum and angular momentum vectors. As Eq ([5]) is of 
this form, four of the integrals of motion are already known and arise from separability. The 
puzzle is that there is a fifth integral whose form is unknown and which does not arise from 
separability. 

In the case where one of the barriers, say that in the x = y = plane, vanishes, the fifth 
integral is known to be [f| 

k k-i k 2 
h — + — 

2r x 2 y 2 

and follows from separability in the rotational parabolic coordinate system. If the fifth 
integral for the general problem (pjj) is quartic, it must reduce to the above integral in the 
limit k 3 — > 0. This suggests taking the ansatz 

h = (l iP2 - Pl L 2 - 2z (- A + h 2 + ^ ^ + hg{ x, p) (14) 

where g(x, p) is a function of both position and momentum yet to be determined. Requiring 
the Poisson bracket of I4 with the Hamiltonian to vanish leads to a solution for g. So, we 
arrive at an isolating integral of the form 

( 1 \ ( k k\ ko k^\\ 2 Ik-i ,0 . , 

^((Lxp),-^--^^^)) +^(-P) 2 (15) 

It is not yet proven that this is a functionally independent integral, as it is possible to 
construct an infinite number of quartic integrals from combinations of the four existing 
quadratic integrals. To test for functional independence, the 5x6 Jacobian 

d(E,h,I 2 ,I 3 ,h) 
d{x hPi ) 

can be constructed and shown to be of rank 5. Thus, the integral given in Eq (fl5l) is the 
fifth functionally independent isolating integral of motion for the Hamiltonian. 



In the case when the ki are all zero, this integral reduces to the z-component of the 
Laplace- Runge-Lenz vector. In fact, if we work through the same derivation but make the 
cyclic permutations x — > y — > z, two more integrals are obtained, They are the equivalents 
of Eq ffTBT) with the coordinates permuted, and reduce to the x— and y— components of the 
Laplace- Runge-Lenz vector. As is expected they are not functionally independent, and this 
is easily demonstrated by including them in the Jacobian and noting that it remains of rank 
5. The three components are related through 

h x + h y + K = ±E{h + h + h + h + k 2 + k 3 ) + k 2 (17) 

where I± z is the integral given by Eq (1151) and and the cyclicly permuted versions. 



IV. ACTION- ANGLE VARIABLES 



It is also possible to solve this general problem in action-angle variables. This is worth- 
while as it provides insight into the relationship our of our new superintegrable Hamiltonian 
with the Kepler problem. Following [14J], the actions can be shown to be 



and hence 




As expected, the Hamiltonian depends on the actions only through the combination J r + Jg + 
Jfa implying that the three frequencies of the classical motion are the same. If (w; r , u^, u^) 
are the angles conjugate to (fl9|) . then we can make a canonical transformation to new 
action-angle coordinates (Ji, J 2 , J3,Wi,w 2 ,w 3 ), using the generating function 

F = {wtj, - w e )Ji + (w e - w r )J 2 + w r J 3 (20) 

The new actions are related to the old via 

Ji — J<t>i J2 — Je + J 4,1 J3 — Jr + Jo + J(j)i (21) 
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and so the Hamiltonian becomes 

-?k 2 TT 2 

H = E= — (22) 

and depends on only one of the new actions. Using Hamilton's equations, we see that the 
angle W3 increases linearly with time, whilst the angles w% and w 2 are the additional integrals 
of motion. They can be found by explicit construction of Hamilton's characteristic function 
S 

S = p<j,d(j) + / p e d9 + / p r dr (23) 



followed by use of the equations 10, = dS/ dJi (see [14( for the equivalent calculation for the 
Keplerian potential). We find that 

1 (A 2 cos 2(f) - ki + k 2 

Wi = — arccos — 

4vr ^(A 2 - fci + fc 2 ) 2 - 4A 2 fc 2 

1 / 2A 2 cot 2 6 - B 2 + A 2 + k 2 . 

arccos — . (24) 

4vr \V(B 2 - A 2 - k 2 ) 2 - AA 2 k 2 ' 

1 f B 2 cos 26 + A 2 - k 

w 2 = — arcsm 



where 



4tt \^(A 2 + B 2 -h) 2 -AA 2 B< 

1 . / kr-2B 2 \ 

— arcsm — = (25) 

2tt \rVk 2 + AB 2 Ej 

A = + V / ^+ Vh = VT2 (26) 

2y 2tt 

B = + Vh + ^fk 2 + v^= s/h (27) 

2y 2ir 

Note that in the case that k\ = k 2 = k^ = (the Kepler problem), it is usual to introduce 
the inclination i = arccos(A/i?) of the orbital plane. Tnen, the angle W\ reduces to 

w\ = — (<f) — arcsin (cot 6 cot i)) (28) 
2tt 

which is the longitude of the ascending node. The second angle w 2 is easiest evaluated in 
the orbital plane with polar coordinates (r, <p) and becomes 

W 2 = 7T (if ~ V^lan) - (V 9 - fperi)) = W. (29) 

Zix 

The first integral is therefore the angular difference between the orbital position and the 
longitude of the ascending node, the second the angular difference between the orbital po- 
sition and the periapse. Thus, w 2 reduces to u>, the longitude of the periapse in the Kepler 
problem. 



V. SUMMARY AND CONCLUSIONS 



We have found a new superintegrable Hamiltonian, which is a generalization of the well- 
known Kepler problem. There are five isolating integrals of the motion, namely the energy 
and generalizations of the componenets of the angular momentum and Laplace-Runge-Lenz 
vectors. Intriguingly, some of the integrals of motion are quartic in the momenta and do not 
arise from separability of the Hamilton- Jacobi equation. 

There are three interesting questions which merit further research. First, it is clear that 
the N degrees of freedom Hamiltonian 



2N — 1 functionally 
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16] that the ad- 



is also superintegrable. It would be interesting to find the complete set o 
independent integrals of motion in this case. Second, it is well-known 
ditional integrals of motion in the Kepler problem arise from the existence of the dynamical 
symmetry group SO (4). It would be interesting to understand the group theoretic inter- 
pretation of the integrals of motion discussed in this paper. Third, although the reduction 
technique we used to generate the superintegrable potential in Section [III is simple, it can be 
made to do some more work. For example, it is also clear that the Hamiltonian [c.f., eq (jBJ)] 

H = ~ (pi + vl + Pi) + fx 2 + m 2 y 2 + n 2 z 2 + -± + -§ + -§. 

always has 5 independent integrals of motion as well, whose form remains to be estab- 
lished. Perhaps all superintegrable potentials in three degrees of freedom can be viewed as 
projections of higher dimensional Keplerian or harmonic oscillator motion? 
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